Self-organization of the MinE ring in subcellular Min oscillations 
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We model the self-organization of the MinE ring that is observed during subcellular oscillations of 
the proteins MinD and MinE within the rod-shaped bacterium Escherichia coli. With a steady-state 
approximation, we can study the MinE-ring generically - apart from the other details of the Min 
oscillation. Rebinding of MinE to depolymerizing MinD filament tips controls MinE ring formation 
through a scaled cell shape parameter f . We find two types of E-ring profiles near the filament 
tip: a strong plateau-like E-ring controlled by ID diffusion of MinE along the bacterial length, or 
a weak cusp-like E-ring controlled by 3D diffusion near the filament tip. While the width of a 
strong E-ring depends on f , the occupation fraction of MinE at the MinD filament tip is saturated 
and hence the depolymerization speed do not depend strongly on f. Conversely, for weak E-rings 
both f and the MinE to MinD stoichiometry strongly control the tip occupation and hence the 
depolymerization speed. MinE rings in vivo are close to the threshold between weak and strong, 
and so MinD-filament depolymerization speed should be sensitive to cell shape, stoichiometry, and 
the MinE-rebinding rate. We also find that the transient to MinE-ring formation is quite long in 
the appropriate open geometry for assays of ATPase activity in vitro, explaining the long delays 
of ATPase activity observed for smaller MinE concentrations in those assays without the need to 
invoke cooperative MinE activity. 

PACS numbers: 87.17.Ee, 87.16.A-, 87.16.dr 



I. INTRODUCTION 



The oscillation of the proteins MinD and MinE from 
pole to pole of individual cells of the bacterium Es- 
cherichia coli is used to localize cellular division to mid- 
cell flj. One cycle of the oscillation, lasting approxi- 
mately one minute, starts with ATP-associated MinD 
binding to the bacterial inner membrane and polymer- 
izing into helical filaments 0, i, H (see also (1). This 
occurs at alternating poles of the bacterium, with the 
MinD forming a polar "cap". MinE is recruited to the 
membrane-bound MinD, where it forms a distinctive "E- 
ring" [1, 0, [1] at the edge of the MinD cap by accu- 
mulating near the MinD filament tips [2]. Because the 
rate of hydrolysis and subsequent release of ATP-MinD 
is stimulated by MinE [1, 0, @] , the E-ring drives depoly- 
merization of the MinD filament which allows the oscil- 
lation to proceed. The depolymerization occurs with an 
approximately fixed E-ring width and speed along the 
cell axis 0, indicating an approximate steady-state 
during this part of the Min oscillation. However, little 
is known about the mechanism of E-ring formation, its 
detailed structure, or how important it is for Min os- 
cillations. Indeed, Min oscillations have been observed 
without prominent E-rings 10]. 

Most models proposed for Min oscillation do not have 
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explicit MinD filaments [Ul, [ij, [ij, UJ, ll5|, though 
they do have E-rings. Recently, several models of Min 
oscillations that include exp licit MinD polymerization 
have been proposed [IB, IITL ITsI . [l9l | , two of which dis- 
play strong E-rings that track the tips of depolymeriz- 
ing MinD filament caps with constant speed and width 
[H, [3- In these models, E-rings are the result of MinE 
polymerization either orthogonal to or along 1^, 
MinD filaments. While MinD polymerization has been 
observed in vitro 0, there have been no reports of 
MinE polymerization in the experimental literature. In- 
deed, the faint MinE "zebra-stripes" associated with the 
MinD zones adjacent to the MinE ring 0, d, [l3| seem 
to imply sparse lateral binding of MinE to the body of 
MinD filaments - not MinE polymerization. 

In this paper, with both stochastic 3D simulations, 
and a deterministic ID model, we show that local (non- 
polymeric) rebinding of MinE released from depolymeriz- 
ing MinD filament tips is sufficient for E-ring formation. 
We impose and characterize a dynamical steady-state of 
an E-ring on a depolymerizing semi-infinite MinD fila- 
ment in order to address the approximate steady-state 
speed and width of the E-ring in vivo 0, HI- We in- 
vestigate the roles of spatial dimension, cell length, and 
radius, and of multiple MinD filaments and their heli- 
cal pitch. We estimate the timcscale of E-ring formation 
and obtain results consistent with the significant delay 
before ATPase activity seen with small MinE concentra- 
tions and large MinD membrane coverage in vitro [3, ^. 
Finally, we discuss how competition between the intrin- 
sic and the MinE-stimulated ATPase activity of MinD 
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FIG. 1: Fractional occupancy of MinE on the MinD filament 
p vs. distance along the bacterial axis z for the 3D stochastic 
model (continuous lines) and the ID model (dashed lines) for 
parameters typical of E. coli: L = 2 /xm, R = 0.5 /im, and 
po = 0.35. The MinE binding parameter was era = 
while for the ID model / = 0.06 was used. One MinD fila- 
ment supports either (A) a strong plateau-like E-ring for pitch 
p = 0.45 /im or (B) a weak cusp-like E-ring for pitch p — oo 
(straight filament). The width W of the strong E-ring, given 
by p{W) = (1 -1- po)/2 is indicated. The inset illustrates the 
cylindrical geometry of the 3D model, showing the underlying 
helical MinD filament with its depolymerizing tip at a = 0. 
The helical pitch p is indicated. 



MinD monomers; if not it is removed from the system at 
z = ±L. This open boundary condition represents the 
sinks for MinE provided by other MinD in the system. 
Dcpolymcrized MinD is removed from the system with- 
out further interaction, reflecting the nucleotide exchange 
needed before MinD rebinding is possible. This dramat- 
ically simplifies our model, since we may then explicitly 
consider only MinE dynamics on an implicit MinD fil- 
ament. Both the boundary conditions and the neglect 
of depolymerized MinD will be addressed again in the 
discussion. 

In order to study a steady-state E-ring, we keep the 
filament tips centered at z = - the "tip-frame" . In the 
tip-frame, bound MinE move along MinD filaments at a 
constant depolymerization speed v while new monomers 
of MinD are introduced at z = L decorated with MinE 
with a constant probability po (determined by the rela- 
tive cellular amounts of MinE and MinD particles). [Ef- 
fectively we are studying semi-infinite MinD filaments 
under the approximation of uniform MinE binding for 
z > L.] In the steady-state, the fraction of MinE released 
by the depolymerizing MinD filament tip that reach the 
absorbing boundaries will then be po/pup, where pup 
is the fractional MinE occupation of the filament tip. 
The depolymerization speed v can be determined self- 
consistently by pup, though we will see below that v is 
small and can be practically ignored in terms of the E- 
ring structure. 



controls the instability that leads to the initial formation 
of the E-ring from a uniformly decorated MinD filament. 

Qualitatively, we predict that the width of MinE- 
rings will increase as the MinD-filament depolymeriza- 
tion speed is increased through manipulation of cell 
shape, MinD to MinE stoichiometry, or mutations that 
affect the MinE binding rate to MinD. Eventually, the de- 
polymerization speed will saturate but the E-ring width 
can still grow. Conversely, as the depolymerization speed 
is decreased, MinE-rings will undergo a transition from 
a plateau-like "strong" E-ring to a cusp-like "weak" E- 
ring. To our knowledge, systematic experimental studies 
of the E-ring width have not yet been done. 



II. E-RING MODEL 

As illustrated in the inset of Fig. [H we represent the 
bacterial geometry as a cylinder of radius R and length 
2L. In the right half (0 < z < L), n filaments of MinD are 
placed on the cylinder, each with the same helical pitch 
p but with random (unbundled) helical phases. MinD 
filaments are composed of monomers of length oq, each 
of which can bind one MinE. We depolymerize MinD 
from filament tips at z = 0, and any released MinE dif- 
fuses in the cylinder interior (cytoplasm) with a diffu- 
sion constant D. Released MinE can bind to unoccupied 



A. Stochastic 3D implementation 

The dimensionless parameter ai = v£/D, the frac- 
tional axial distance one MinE advects at speed v while 
it diffuses a distance i, characterizes the importance 
of the depolymerization speed. Even with £ — Ap,m, 
V = Omp,m/s [^,11, and L» = lOpm^/s ai ^ O.OI is 
small and depolymerization is slow compared to diffusion. 
Accordingly, our stochastic 3D model quasiadiabatically 
follows each released MinE until it either rebinds or is 
removed from the system before allowing further depoly- 
merization. Each MinE diffuses by taking a randomly 
oriented step of fixed length 6 every timestep At, where 
D = S^/{6At). Diffusing MinE binds to a free MinD with 
probability PsUck when it hits the bacterial membrane 
within a distance rbind of the MinD. We take rbmd = 0,0 ■ 
This leads to an effective binding rate of u^p^^iocah where 
Pa.iocai is the local bulk concentration of MinE and the 
bulk reaction rate 0-3 = SnDrl-^^PsUck / {"^S) ■ We take 
CT3 — 0.3 pm? /s (this is approximately the threshold be- 
tween strong and weak E-rings given the cell geometry, 
see below) . The steady state reached after successive de- 
polymerization steps is independent of small 5 if we vary 
Pstick with 5 to keep (T3 constant. 
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B. Analytic ID treatment 

We also study a deterministic ID model that exactly 
corresponds to the 3D stochastic model in the limit 
i? <C ao- This enables us to explore the role of spa- 
tial dimension and stochastic effects in the E-ring, and 
also helps us to identify the combinations of parameters 
that control the E-ring structure. Our ID model tracks 
both the linear density of bound MinE (B) and of freely 
diffusing MinE {F): 

B-vB' = aiF{B,nax - B) , for z > (1) 
F-vF' = DF" -aiF{Bmax- B)+vB{0)5{z),{2) 

where the dots and primes indicate time and spatial 
derivatives, respectively. For z < there are no fila- 
ments so i? = 0. For z > 0, the linear density of potential 
binding sites (i.e. of MinD) is B^ax = and the ID 
rebinding rate is cti. The v dependent terms on the left 
side of the equations represent advection of bound MinE 
in the tip frame, while on the right of Eqn. [5] is a source 
term due to MinE release at the depolymerizing filament 
tip. If we rescale all lengths by L (so z = z/L) and define 
dimensionless fields B = B/B^ax and F = DaF/{vL), 
then we can consider the scaled steady-state equations: 

B' = -CTiF(l - S) , for i > (3) 
F" = -ulF' - B' - B{0)S{S). (4) 

The boundary conditions are -F(±l) = and B{1) — po- 
The behavior is controlled by the dimensionless param- 
eters ai = aiL'^ / (Da) and ul — vL/D, as well as by 
pt). We integrate Eqn. [5] for z < where B = 0, and 
impose flux conservation of MinE at the boundaries with 
F'{—1) — F'{1) = Po- For z > the equations are then 
integrated numerically to find the steady-state. 

Following the discussion of the stochastic 3D imple- 
mentation, we expect to be small, and anticipate that 
it is irrelevant for the E-ring structure - leaving only po 
and (Ti as relevant control parameters. Nevertheless, the 
ID treatment allows us to explore this assumption. We 
find that a < 0.05 does not change the observed E-ring 
steady-state structure by eye, while we expect ai « 0.01 
at room temperature in vivo — and even lower values 
for weak E-rings. The four-fold speedup observed for the 
Min oscillation at body temperature [2l| puts the de- 
polymerization speed (i.e. a) closer to, but still under, 
relevance with respect to the structure of the steady-state 
MinE ring. 

III. RESULTS 

We can compare results of our ID deterministic model 
with our 3D stochastic model using F = TrR^ps^av, where 
Ps.av is the bulk density averaged over the bacterial cross- 
section. Then the ID and 3D binding rates of MinE are 
related by ai = a^f /{ttR^), where / = pz.iocail Pz^av We 



expect that / will vary with distance from the filament 
tip due to local release at the tip followed by diffusion 
and capture. We find / < 1 away from the filament tip 
due to rebinding to the MinD filament, and we expect 
/ > 1 at the filament tip due to local release from the 
depolymerizing tip. Effects of multiple filaments (n) and 
filament pitch (p) can be included in the ID model by 
using the MinD monomer spacing projected along the 
bacterial axis o, where 

a = ao/(?^v/l + 47r2i?7p2)_ (5) 

Differences between the two approaches are either due to 
the ID vs. 3D geometry or due to the deterministic vs. 
stochastic nature of the models. 



A. Strong and weak E-rings 

Fig. [1] illustrates the fractional occupation p (equiva- 
lent to B in the ID model) of MinE binding sites on the 
MinD filament vs. distance z along the bacterial axis. 
Occupation monotonically decreases from the tip value, 
ptip = /o(0), due to local rebinding of MinE following 
depolymcrization from the tip. Following the quantifica- 
tion of Shih e< al there are a few thousand MinD 
monomers within a typical bacteria. With L = 2/im 
and ao = 5nm [3^, '3| , they can be arranged either in one 
single helical filament (with p « 0.45/im 0) or about 
7 straight filaments (with p = 00). In either case, we 
find (A) a "strong" E-ring [n — \ shown) with pup ~ 1 
and a plateau shape of the density profile near the tip. 
With only one straight filament (B), we find a "weak" E- 
ring with enhanced density at the tip but no saturation 
(ptip < 1) and no plateau. Strong or weak E-rings have, 
respectively, negative (p"(0) < 0) or positive (p"(0) > 0) 
curvature at the tip. 

The ID model profiles reasonably match the 3D re- 
sults away from the filament tips, using / = 0.06. This 
best value of / depends on rjyind- Using the same / 
near the tips, the ID model systematically underesti- 
mates the fractional occupation. This implies that a 
larger / = pz^iocai/ P'i,av is appropriate there, in agree- 
ment with the increased likelihood that MinE will be 
found near the tip shortly after it is released at the tip. 

B. Scaling collapse of E-ring width 

For both strong and weak E-rings, we can define the 
width W of the E-ring such that p{W) = {l + po)/2. Mo- 
tivated by the importance of the scaled MinE rebinding 
rate cfi in the ID deterministic equations and by the cor- 
respondence of (Ti and (T3 , we investigated the influence of 
the scaled aspect ratio f = \J~[fd\ = Rj ixDaj oj, on 
the profile shape, as characterized by W jL and by ptip^ 
in Fig.[2]for both the 3D stochastic model (symbols) and 
the ID deterministic model (lines). Two regimes are de- 
marcated by a vertical dashed line: for small f we have a 
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FIG. 2: MinE profile, characterized by W/L and pup, as a 
function of the scaled aspect ratio f. (a, b): Straight filaments 
{p — oo, po — 0.35). ID model results are shown with solid 
lines (using / — 0.06); 3D stochastic results are indicated by 
symbols. Single filaments {n — 1, green data) for L — 1 pm 
(□), L = 2 ^m (o) and L = 3 ^m (A); multiple filaments 
(n —2, or 5, blue data) for L = l^m (v). (c, d): Helical 
filaments (n = 1, po = 0.35) with L/p =20 (□, red), 10 
(o, green), 4 (A, blue) and (o, pink). For all these data, 
(73 — 0.3/im^/s, and R is varied to explore f. Similar results 
are obtained when a 3 is varied. 




r r 



FIG. 3: (a, b): MinE profile, as characterized by W/L and 
ptip obtained by the 3D model (points) and the ID model 
(lines, using / = 0.06), as a function of r for different values of 
stoichiometry; po — 0.2 (□, continuous lines), 0.35 (o, dashed 
lines), 0.50 (A, dotted lines). For each stoichiometry, the 
same collapse as Fig. [2ja,b) is obtained: data are compiled for 
n=l, 2, 3, 4 and 5, L—1, 2 and 3pm, p = 00, (T3 = O.Spm^ / s, 
and R varies to explore f. Similar results are obtained when 
(73 is varied. 



number of MinE and MinD particles) also controls the 
scaling curves of W/L or pup vs. r. Agreement between 
ID and 3D models for small r and scaling collapse are 
preserved for each pQ. 



C. Correspondence with in vivo Min oscillations 



strong E-ring with p"(0) < 0, a saturated tip {pup « 1), 
and good agreement between the ID and 3D models for 
the E-ring width; for larger r we have a weak E-ring with 
p"(0) > 0, Ptip no longer saturated, and a smaller width 
W. 

The agreement between the 3D and ID results for pup 
and W at small f shows that the essential physics of 
strong E-rings is one-dimensional. For small enough R 
the bacterial cross-section is well explored by MinE by 
the time it has diffused to free binding sites a distance W 
from the filament tip. However, by effectively averaging 
the radial profile the ID model systematically underesti- 
mates the occupation fraction near the tip, as seen with 
Ptip in Fig.[2]and also in the profiles shown in Fig.[T] The 
disagreement becomes stronger as f increases, refiecting 
the increasingly 3D character of the stochastic system at 
larger aspect ratios. However, the system still exhibit a 
remarkable collapse for all values of f. This shows that al- 
though the ID model misses important details about the 
tip enhancement, the scaling behaviour of the 3D system 
with straight filaments is similar to the ID model. 

As shown in Fig. ^c, d), f also captures the effects of 
helical MinD filaments. Smaller pitches lead to stronger 
E-rings. However, the 3D stochastic results do not show 
scaling collapse with respect to f as the monomer spacing 
along the filament qq is a relevant length-scale in addition 
to the projected axial monomer spacing a. Since the ID 
model only uses the effective a, it incorrectly exhibits 
perfect scaling collapse. 

As shown in Fig. ^a) and (b), po (the ratio of the 



Experimentally, W/L « 0.3 is observed in rod-shaped 
cells [1, 0, Q I where we take L as half the bacterial length. 
Using Po ~ 0.35, which is consistent with the ratio of 
MinE to MinD if we assume MinE are always dimerized 
jlQ] . then from Fig.[3l[a) we see that W/L « 0.3 is recov- 
ered for f 0.07 — which corresponds to (73 « 0.3pm'^/s. 
(This as is of the same order of ma gni tude as used in 
a number of previous models in 3D [l3l . [l6| and in ID 
[H, 113) [Hi if we assume R = 0.5pm.) Interestingly, 
this indicates that the E-ring of the normal wild-type 
(WT) Min oscillations is a strong E-ring (with a plateau 
of MinE occupation near the MinD filament tip) but near 
the margin between weak (with ptip < 1) and strong. 
This implies (see Eqn. [6] below) that the tip occupation 
Ptip, and hence the depolymerization speed and the oscil- 
lation period, will strongly depend on the stoichiometry 
of MinE to MinD. Since ks/kj ^ 1, changes to pttp even 
at the percent level should be significant. Indeed, MinD 
overexpression leads to a 2.5-fold increase in the period 
[2^ . This also implies from Fig. [Sj^a) that the width of 
the E-ring will strongly depend on the stoichiometry — 
though this has not (yet) been explored experimentally. 
At a fixed stoichiometry of MinE to MinD (po), we ex- 
pect that overexpression of Min will increase the number 
of filaments and/or decrease the pitch. As a result, we 
expect a slightly stronger E-ring, and a slightly faster 
period - as seen [2^ . 

Optically reconstructed E-rings [2] show a plateau-like 
decoration along the MinD filament, consistent with a 
strong E-ring. E-rings have also been seen in long fila- 
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mentous cells [1,0, Hi and exhibit approximately the same 
width W, though both the spacing between MinD caps 
and the cell length are considerably longer in filamen- 
tous cells than in rod-shaped cells. This indicates that 
the effective L may not be determined by cell shape, but 
rather by other processes preserved between rod-shaped 
and filamentous bacteria such as the length of the MinD 
filaments or spontaneous lateral release (without MinD 
hydrolysis) of MinE away from the tip of the MinD fila- 
ment. 

Shih et al. [10] identified MinE point-mutants 
(MinE^^'S'^ and MinE^''^^) that led to fainter E-rings, 
and double mutants (MinE^'^^^/^'^^'^) that resulted in 
most of the MinE being cytoplasmic with no strong fi- 
rings. Assembly and disassembly of MinD polar zones 
continued with no more than doubled periods [l0| — 
too rapid to be explained by intrinsic depolymerization 
alone (in contrast, see [l^). From Eqn. [H] (below) the 
observed disassembly rates would only require a moder- 
ately enhanced pup ~ 0.9, i.e. a weak E-ring. Indeed, 
in all of these constructs there appears to be enhanced 
CO- localization of MinE with the MinD polar zones ■ 
We believe that the lack of visible E-rings in these mu- 
tants can be explained with decreased as (as suggested 
previously by [13|) and/or enhanced spontaneous MinE 
unbinding away from filament tips. Local rebinding of 
MinE near filament tips would still lead to an enhanced 
Ptip- We predict that the oscillation period in these mu- 
tants should be strongly susceptible to the MinE to MinD 
stoichiometry. 



IV. TRANSIENTS 

We may use our models to check that the transients 
before steady-state are fast enough in the context of the 
normal Min oscillation. If we initially decorate the MinD 
filament with MinE monomers released from z = —L 
consistent with MinE released from a different depoly- 
merizing MinD cap, we find (data not shown) an ini- 
tial decoration pattern that has a plateau-like strong fi- 
ring from the beginning (as previously noted [ll]), so 
that we expect rapid fi-ring formation without appre- 
ciable delay during Min oscillations (as also observed 
experimentally 0, Q)- 

A. Transients before the steady-state in vitro 

While delays are not observed for E-ring formation 
during Min oscillations in vivo, significant delays are ob- 
served in vitro. MinD binds to phospholipid vesicles in 
the presence of ATP and undergoes self-assembly, con- 
stricting the vesicles into tubes with diameters on the 
order of 100 nm Q. Electron-microscopy revealed that 
MinD assembles into a tightly wound helix on the sur- 
face of these tubulated vesicles with a pitch (helical re- 
peat distance) of only 5 nm. Hu et al. Q report a sig- 



nificant delay (several minutes) for stimulated ATPase 
activity when small concentrations of MinE were added, 
while this delay vanished for larger MinE concentrations. 
Similar delays were seen in vitro by Suefuji et al. 
Furthermore, the eventual steady-state ATPase activity 
was smaller for smaller concentrations of MinE 0, |J. 
This has led to the hypothesis of explicit cooperativ- 
ity of MinE binding, which has then been explicitly in- 
cluded in reaction-diffusion models [l3, and in MinE 
polymerization in models with MinD polymers [H, [l^ . 
Here we show that our stochastic model for the MinE 
ring, with no explicit MinE cooperativity, can recover 
the MinE concentration dependent ATPase delays and 
activities observed in vivo. We conclude that MinE co- 
operativity is not needed to explain the in vitro results, 
apart from cooperative effects that arise implicitly from 
the self-organization of the MinE ring. 

We use an "inside-out" open geometry corresponding 
to what is reported in vitro 0], with a narrow phospho- 
lipid cylinder that is tightly wound by MinD filaments. 
MinE, when released by a depolymerizing filament tip, 
will diffuse outside the cylinder. We consider a heli- 
cal MinD filament of radius R = 50 nm and pitch 5 
nm (equal to oq). Upon MinD depolymerization, we al- 
low any released MinE to diffuse until either it binds to 
an available MinD binding site or it is absorbed by the 
boundaries at z — ±L. We impose reflecting boundary 
conditions at r = i?, but otherwise allow MinE to diffuse 
freely for r > R. Our stochastic 3D model is otherwise 
the same as before though with an emphasis on the tran- 
sients approaching steady-state. 

The transient to steady-state is shown in Fig. Ufa), 
with the fractional occupation of MinE at the MinD fil- 
ament tip {ptip) shown as a function of the number of 
depolymcrized monomers from the filament tip, TV. The 
MinE occupation fraction at the MinD filament tip is ex- 
perimentally observable through the ATPase activity (i.e. 
the MinD depolymerization rate). The initial condition 
is a uniform occupation po, corresponding to an initially 
random binding of MinE on the MinD filament. The 
larger po is, the shorter the transient and the stronger 
the eventual steady-state pup. Significant enhancement 
of Ptip is obtained even for small fractions of MinE. For 
Pa It 0-2 we see that pup > 0.8, though, as shown in 
Fig. H^b) , strong E-rings are predicted only for very large 
stoichiometry (po ^ 0.8). The inside-out in vitro geome- 
try includes some small radius features (the helical wind- 
ing of the MinD filament) and some large radius features 
(no closed boundary at large r). The tight helical wind- 
ing of the MinD filament contributes to long transients, 
while the semi-infinite radial geometry contributes to the 
weak E-ring for small and moderate po- 

To convert the number of depolymerization steps N to 
a time t{N) we need to sum the average time for each 

step, which will depend on pup'. t{N) — J2n=i ^^i"^) 
where, 

At(ri) = ptip{n)/ks + (1 - Ptip{n))/ki. (6) 
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FIG. 4: Transients and E-ring structure for an "inside-out" 
open geometry appropriate for m vitro experiments, where 
a MinD filament is tightly wound on the outside of a cyhn- 
der of small radius [R — 50 nm) with open boundaries at 
R = oo. (a) Evolution of pup as a function of the number of 
depolymerization steps A'^ (measured in thousands) after the 
uniform intial conditions for po equal to 0.8 (solid, red), 0.4 
(long dash, green), 0.3 (short dash, blue), and 0.2 (dotted, 
pink); (b) steady-state p{z) as a function of axial distance z 
along the helical axis for the same po- 
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FIG. 5: For the same inside-out in vitro geometry described 
in the previous figure. Cumulative ATP-ase activity N{t) 
(measured in thousands of depolymerization steps ) versus 
time t for various po- Asymptotic behavior are plotted as 
thin dotted lines. 



The timesteps are determined by kj when the tip of the 
MinD polymer is unoccupied by MinE and fcg when it is 
occupied. Using ks/kj ~ 20^ and ks — l/(20ms) given 
by the maximal depolymerization speed in vivo (assum- 
ing Ptip « 1, with a strong E-ring) [2q . we plot the cumu- 
lative total ATPase activity N{t) (equal to the number 
of depolymerization steps) vs. elapsed time t in Fig. [51 



The stoichiometric ratio of MinE to MinD corresponds 
to Po if the MinE mostly binds to available MinD be- 
fore depolymerization proceeds significantly. For small 
amount of MinE (typically po < 0.3) we obtain a sig- 
nificant delay of about 5 minutes, corresponding to the 
ATPase delay seen in vitro 3 ; a-nd for larger MinE 
amounts (po going to 1) the delays decrease towards zero 
also in agreement with in vitro studies. When the steady- 
state Ptip is reached, the ATPase rate will also be in 
a steady-state as indicated by the linear asymptotes in 
Fig. O Since ptip can be large even for smaller po, we 
expect the ATPase rates to be comparable for moderate 
or larger po, as seen in vitro [1, 01 • For smaller po the 
steady-state ATPase activity is reduced, as also observed. 

We conclude that the delay of ATPase activity seen 
in vitro is determined by the time needed to reach the 
steady-state ptip- We see that it is considerably longer 
in an open than in a closed geometry. Our local re- 
binding model recovers the delays seen in vitro with- 
out any exp licit MinE cooperativity (see, conversely, 

iiHMEHl). 



B. E-ring instability 

In the tip-frame, the MinD filament tip is bistable dur- 
ing Min oscillations ^] and the formation of the E-ring 
switches the filament tip between polymerization and de- 
polymerization. While long transients for this switch- 
ing are not expected during Min oscillations in vivo be- 
cause of initially non-uniform tip decoration [23j , we may 
ask about the transient to form the E-ring from a non- 
oscillating state — such as seen experimentally after ex- 
posure to high levels of extracellular cations We 
consider a MinD filament that is initially uniformly dec- 
orated with MinE. To tractably include the MinD poly- 
merization dynamics, we use a uniform (mean- field) bulk 
MinD density pjj. Because we are interested in the ini- 
tial slow stages of E-ring formation, we consider MinE 
binding only near the tip with occupation fraction ptip 
(initially equal to po) 

The net polymerization rate of a MinD filament is i? = 
k+PD — {ptipks + i^ — Ptip)ki), where po is the bulk MinD 
monomer concentration and A:+ controls MinD monomer 
addition. Depolymerization of n monomers from a single 
tip will enhance ptip due to local rebinding of MinE, so 
that dR/dn — k+/V — [ks — ki)dptip / dn for cell volume 
V . The depolymerization time per monomer is Ai « 
1/ki for an initially weak E-ring (with ptip small), and 
the change in tip occupation in one depolymerization step 
will be proportional to both the number of MinE released 
{ptip) and the locally available binding sites (1 — ptip), so 
that 

dR 

— = k+ki/V - A{ks - ki)kiptip{l ~ pttp), (7) 

where the constant A is the fraction of MinE that re- 
bind to available sites at the filament tip. For kg suf- 
ficiently greater than kj this represents an instability 



7 



[dR/dt growing more negative with time) that will lead 
to E-ring formation. We therefore expect that both a sig- 
nificant difference between intrinsic and stimulated AT- 
Pase activity of MinD and significant intrinsic ATPasc 
activity are needed for E-ring formation, and hence for 
the initiation of Min oscillations. 

We have neglected any lateral unbinding of MinE from 
the MinD filament, which will kill the instability if dR/dt 
is small enough. We also neglect the presence of other 
MinD filament tips, which will buffer the bulk MinD den- 
sity and reduce the effect of the fc+ term in Eqn. [71 These 
effects will shift the threshold, but will not change the 
presence of the E-ring instability. 

Since pup — po initially, we also predict from Eqn. [7| 
that both low and high proportions of MinE to MinD will 
also preclude Min oscillations by making the MinD fila- 
ment tip initially stable against depolymerization. How- 
ever, using k+ = 100/{fj.Ms) [l^, A « 1, and V = l^m^ 
we estimate a tiny stoichiometry threshold of 0.003 (for 
Po or 1 — Po ) • While our predicted stoichiometry thresh- 
olds are unlikely to be relevant in vivo, they may be 
approachable in vitro. We also note that initially slow 
E-ring formation dynamics near the instability threshold 
should be observable when Min oscillations are restarted 
after being halted ^6| . 

Previous models of the full Min oscillation have found 
limiting MinE: MinD stoichiometrieSjeither both low and 
high [li, [H [li] or just high [H, Sufficiently low 

stoichiometrics may not have been explored in the later 
models. Conversely, Min oscillations have always been 
seen in vivo with moderate stoichiometry changes [2^ . 
It would be desirable for a more systematic exploration 
of the role of stoichiometry on Min oscillations, given 
the predicted stoichiometry limits for the existence of 
oscillations predicted in this and other models. 



V. DISCUSSION 

We have presented a model of the self-assembly of the 
MinE-ring within single E. coli bacteria, without invok- 
ing either MinE cooperativity or MinE polymerization. 
We highlight the difference between strong E-rings, with 
Ptip ~ 1, essentially ID physics and a maximal depoly- 
merization speed, and weak E-rings with p^p < 1 that 
have 3D physics with depolymerization speeds that sen- 
sitively depend on the parameters, and especially on 
the amount of MinE in the cell. In contrast to pre- 
vious filamentous models that had only strong E-rings 
[isl . [igj . our model shows how changing the stoichiom- 
etry of MinE and MinD can change the oscillation pe- 
riod through the depolymerization speed of MinD fila- 
ments. MinE-rings in non-polymeric reaction-diffusion 
models [ll|, [l^, [IE 113 follow essentially our local 
rebinding mechanism in the ID regime, but will devi- 
ate from polymeric models for weaker E-rings in the 3D 
regime where the monomer scale Oo enters. Since the 
experimentally measured E-ring width indicates that E- 



rings in vivo are close to the threshold between weak 
and strong, the detailed response of the E-ring structure 
(i.e. the width W, or the depolymerization speed via the 
tip occupation pup) to experimental manipulations that 
change the oscillation period (stoichiometry through po 
or, e.g., [1^) is unlikely to be correctly captured by ID 
or non- filamentous models. 

We have explained the anomalous delays of MinE stim- 
ulated MinD ATPase activity seen in vitro [1, 01, and 
have also identified an instability of MinE ring forma- 
tion that is required to develop from a disordered initial 
state to the full Min oscillation. We have shown that 
MinE-ring structure and dynamics can be treated inde- 
pendently of a full Min oscillation model. The instabil- 
ity to E-ring formation, and subsequent MinD filament 
depolymerization, that we identify neither depends on 
nor determines the spatial pattern of Min oscillation - 
which could be selected by either diffusion and rebinding 
of MinD or by phospholipid heterogeneities [2^ . 

We have constructed our E-ring model to obtain a 
steady-state. The steady-state is formed by balancing 
the MinE entering the system as a bound fraction po on 
the MinD filament with the MinD lost by diffusing across 
the open boundaries at z ± L. Other geometries, such as 
an open boundary at z = —L and closed at z = L, or a fil- 
ament tip placed asymmetrically (away from z = 0), will 
also lead to a steady-state E-ring that should be qualita- 
tively similar to the one we have described. An extreme 
example of this is the inside-out geometry we used to 
describe in vitro ATPase experiments. What we have ac- 
complished is to characterize the steady-state, and use it 
to explore the effects of cell-shape, helical pitch, MinE 
rebinding rate, and stoichiometry on the E-ring struc- 
ture. Our model is expected to be a generic part of full 
oscillation models that exhibit E-rings. 

It is worth speculating on how our simplified E-ring 
model would be modified by possible additional ingredi- 
ents within a full model of the Min oscillation. (1) We do 
not expect that filament cutting (see e.g. [l^, [l3]) will 
qualitatively affect our results, though it would lead to 
many more free ends and faster depolymerization. The 
MinE ring would still only be expected to form near the 
very end of the MinD filament, and significant depoly- 
merization would only occur within its width W from 
the end. Similarly, our results should apply to models 
without filaments (see e.g. [H [H, Q [11, El ) . In that 
case, we expect that our analytic ID treatment to be 
a better approximation due to the absence of an intrin- 
sic monomer spacing oq that is relevant near the fila- 
ment tip. (2) We expect that lateral release of bound 
MinE away from filament tips, without associated cut- 
ting, would affect the E-ring profile a distance £ — vt 
away from the tip (where v is the depolymerization rate, 
and is the lateral release rate). This can be crudely 
included in our model by placing our boundary condi- 
tions at i w £. (3) We have neglected the rebinding of 
MinD to the filament tip. We would expect rebinding to 
"poison" the E-ring by significantly reducing the depoly- 
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merization rate - which would allow further rebinding. 
This appears to be observed in the occasional E-ring re- 
versal in vivo [3, [13] ■ While interesting, poisoning ap- 
pears to be typically avoided during Min oscillations — 
perhaps by filament cutting or by lateral MinE release 
and re-binding, neither of which have been experimen- 
tally characterized — and so we are justified in neglect- 
ing it for steady-state E-rings. Poisoning may however 
weaken the E-ring instability described by Eqn. [3 and 
this deserves further study. The next step is to develop 
a full 3D Min oscillation model with MinD filaments but 
without MinE polymerization. 

Previous work has considered the steady-states of semi- 
infinite filaments with tip-directed depolymerization en- 
hanced by bound motors (in this paper, bound MinE) 
[2^. That work used a uniform (mean-field) cytoplas- 
mic motor distribution, and obtained tip-enhanced motor 
density by a combination of diffusion and directed motion 
along the filament together with a "processivity" reten- 
tion probability p for motors at the depolymerizing tip. 
In contrast, in our model MinE remains immobile on the 
filament. [Note that advection (v) represents the drag- 
ging of MinE along with the MinD filament, not motion 
with respect to the filament.] Furthermore, we explicitly 



consider the cytoplasmic MinE random-walk or diffusion 
upon release from the filament tip. While this does lead 
to implicit processivity (local retention of MinE), it also 
correctly allows for rebinding of MinE away from the fil- 
ament tip. This physical modeling of the cytoplasmic 
MinE allows us to consider, e.g., the 3D vs. ID cross- 
over, realistic transients for the inside-out in vitro geom- 
etry, and the E-ring width. Note that the enhanced local 
rebinding of MinE to the MinD filament upon release 
is related to ligand rebinding (see, e.g., [s^]), and sim- 
ilar dimension and geometry dependent effects are seen 
there. 
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